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Abstract 

In this paper we apply Clark-Ocone formula to deduce an explicit 
integral representation for the renormalized self-intersection local time 
of the d-dimensional fractional Brownian motion with Hurst param- 
eter H € (0, 1). As a consequence, we derive the existence of some 
exponential moments for this random variable. 

1 Introduction 

The purpose of this paper is to apply Clark-Ocone's formula to the renor- 
malized self-intersection local time of the d-dimensional fractional Brownian 
motion. As a consequence, we derive the existence of some exponential 
moments for this local time. 

A well-known result in Ito's stochastic calculus asserts that any square 
integrable random variable in the filtration generated by a (i-dimensional 
Brownian motion W = {Wt,t > 0} can be expressed as the sum of its ex- 
pectation plus the stochastic integral of a square integrable adapted process: 



" poo 

F = E{F) + Y J u\t)dWt. 
7Ti Jo 



The process u is determined by F, except on sets of measure zero. In this 
context, Clark-Ocone formula provides an explicit representation of u in 
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terms of the derivative operator in the sense of Malliavin calculus. More 
precisely, if F belongs to the Sobolev space ID 1 ' 2 , then u L (t) = E(D\F\Tt)-, 
where D l denotes the derivative with respect to the ith component of the 
Brownian motion and {J~t, t > 0} is the filtration generated by the Brownian 
motion. Extensions of this formula have been developed by Ustiinel in |17j . 
and by Karatzas, Ocone and Li in |12j . Clark-Ocone formula has proved to 
be a useful tool in finding hedging portfolios in mathematical finance (see, 
for instance, [TT]). 

The fractional Brownian motion on M. d with Hurst parameter H G (0, 1) 
is a d-dimensional Gaussian process B H = {B^,t > 0} with zero mean and 
covariance function given by 

E(B^B^) = 5 -f{t 2H + s 2H -\t- s\ 2H ), (1.1) 
where i,j = l,...,d,s,t>0, and 

r i if i = j 

is the Kronecker symbol. Assume d > 2. The self-intersection local time of 
B H is formally defined as 

L= f T !\{B? -Bf)ds, 
Jo Jo 

where <5o is the Dirac delta function. It measures the amount of time that 
the process spends intersecting itself on the time interval [0,T]. Rigorously, 
L is defined as the limit in L 2 , if it exists, of L e = f Q J Q t p £ (Bl 1 — B^)dsdt, 
as e tends to zero, where p e denotes the heat kernel. 

For H = |, the process B is a classical Brownian motion and its self- 
intersection local time has been studied by many authors (see Albeverio et 
al. p], Calais and Yor [3], He et al. [6], Hu [7j, Imkeller et al. [TO], Varadhan 
|18j . Yor [20], and the references therein). In this case, if d = 2, Varadhan 
[18] has proved that L £ does not converge in L 2 , but it can be renormalized 
so that L e — E(L £ ) converges in L 2 as e tends to zero to a random variable 
that we denote by L. This result has been extended by Rosen [TO] to the 
case He (h, |) (still when d = 2), and by Hu and Nualart in [9], where 
they have obtained the following complete result on the existence of the 
self-intersection local time of the fractional Brownian motion: 

(i) The self- intersection local time L exists if and only if Hd < 1. 
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(ii) If Hd > 1, the renormalized self- intersection local time L exists if and 

3 
2' 



only if Hd < 3 



An important question is the existence of moments and exponential mo- 
ments for the (renormalized) self-intersection local time. Along this direc- 
tion, Le Gall [13] proved that for the planar Brownian motion, there is 

a critical exponent Aq, such that -ElexpAL) < oo for all A < Aq, and 



E ^exp XLj = cxd if A > Aq. Using the theory of large deviations, Bass and 

Chen proved in [2] that the critical exponent Ao coincides with A~ 4 , where 
A is the best constant in the Gagliardo-Nirenberg inequality. 

Clark-Ocone formula seems to be a suitable tool to analyze the renormal- 
ized self-intersection local time, because in this formula we do not take into 
account the expectation of the random variable. The fractional Brownian 
motion can be expressed as the stochastic integral 



Bf = [ t K H (t,s)dW s 
Jo 



of a square integrable kernel Kn(t, s) with respect to an underlying Brown- 
ian motion W. In this way the renormalized self-intersection local time L is 
a functional of the Brownian motion W, and we can obtain an explicit in- 
tegral representation L, in the general case Hd < |. This formula allows us 
to obtain some exponential moments for the renormalized self-intersection 
local time, using the method of moments. 

The paper is organized as follows. In Section 2 we present some prelimi- 
naries on Malliavin calculus and Clark-Ocone formula. Section 3 is devoted 
to derive estimates for the moments of the self-intersection local time in 
the case of a general ci-dimensional Gaussian process, using the method of 
moments. In the case of the fractional Brownian motion, this provides the 
existence of exponential moments in the case Hd < 1. Section 4 contains 
the main result, which is the integral representation of the renormalized 
self-intersection local time of the fractional Brownian motion in the case 
H < min ^t|j, ^py^ . As an application we show that E ^exp L ^ < oo 

if p < ^ [(^ + H) (I — gjf)] . A crucial tool is the local nondeterminism 
property introduced by Berman in [3] and developed by many authors (see 
Xiao [19] and the references therein). 
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2 Preliminaries on Malliavin calculus and 
Clark-Ocone formula 



We need some preliminaries on the Malliavin calculus for the d-dimensional 
Brownian motion W = {Wt,t > 0}. We refer to Malliavin [H] and Nualart 
[15j for a more detailed presentation of this theory. 

We assume that W is defined in a complete probability space (O, J 7 , P), 
and the cr-field T is generated by W. Let us denote by H the Hilbert space 
L 2 (R + ; M. d ) , and for any function h E H we set 



Let S be the class of smooth and cylindrical random variables of the form 



where n > 1, h\, . . . , h n E H, and / is an infinitely differentiable function 
such that together with all its partial derivatives has at most polynomial 
growth order. The derivative operator of the random variable F is defined 
as 



where i = 1, . . . , d and t > 0. In this way, we interpret DF as a random 
variable with values in the Hilbert space H. The derivative is a closable 
operator on L 2 (£l) with values in L 2 (£l;H). We denote by D 1 ' 2 the Hilbert 
spaced defined as the completion of S with respect to the scalar product 



The divergence operator 5 is the adjoint of the derivative operator D. The 
operator 5 is an unbounded operator from L 2 (£l;H) into L 2 (il), and is 
determined by the duality relationship 



for any u in the domain of 5, and F in D 1 ' 2 . Gaveau and Trauber [5] proved 
that 8 is an extension of the classical ltd integral in the sense that any d- 
dimensional square integrable adapted process belongs to the domain of 5, 




F = f(W(h 1 ),...,W(h n )) 



D l t F = Y J ^{W{h 1 ),... 1 W{h n ))h){t) 

j=i J 




E(S(u)F)=E((u,DF) H ), 
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and S(u) coincides with the Ito integral of u: 



d poo 

5(u) = J2 / u\t)dWl 
i=i Jo 



It is well-known that any random variable F G L 2 (Q), possesses a stochastic 
integral representation of the form 

d poo 

F = E(F) + Y, ui (t) dW l 
i=i Jo 

for some d-dimensional square integrable adapted process u. Clark-Ocone 
formula says that if F G B 1,2 , then 

d poo 

F = E(F) + V / E{D\F\Tt)dWl (2.1) 
i=i Jo 



3 Exponential integrability of the self-intersection 
local time 

Suppose that W = {W t ,t > 0} is a ti-dimensional standard Brownian 
motion, defined in a complete probability space (f2, J 7 , P). Suppose that T 
is generated by W. We denote by {F t , t > 0} the filtration generated by W 
and the sets of probability zero. Consider a d-dimensional Gaussian processs 
of the form ^ 

Bt= [ K(t,s)dW s , (3.1) 
Jo 

where K(t, s) is a measurable kernel satisfying J * K(t, s) 2 ds < oo for all 
t > 0. We will assume that K(t, s) = if s > t. 

Fix a time interval [0, T]. We will make use of the following property on 
the kernel K(t, s): 

(HI) For any s,t G [0, T], s < t we have 

t 

K{t,QfdB>ki{t- s) 2H (3.2) 

for some constants k\ > 0, and H G (0, 1). 

Notice that Var (B l t \Fs) = J s * Kit, 9) 2 d6, so condition (HI) is equivalent 
to say that Var (Bi\J- s ^ > k±(t — s) 2H , for each component i = 1, . . . , d. This 
property is satisfied, for instance, in the following two examples: 
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Example 1 Suppose that K(t, s) = (t — s) 2 . Then, we have equality in 
([QD with fei = 

Example 2 Condition (HI) is satisfied by the kernel of the fractional 
Brownian motion, as a consequence of the local nondeterminism property 
(see flUD below). 

We will denote by C a generic constant depending on T, the dimension 
d, and the constants appearing in the hypothesis such as H and k±. 

The self-intersection local time of the process B in the time interval 
[0,T], denoted by L, is defined as the limit in I? as e tends to zero of 



The next theorem asserts that L exists if Hd < 1, and it has exponential 
moments of order -m- 

Theorem 1 Suppose that Hd < 1. Then, the self-intersection local time L 
exists as the limit in 1? of L £ , as e tends to zero, and for all integers n > 1 
we have 




(3.3) 



where p £ denotes the heat kernel 




E(L n ) < C n (n!) 
for some constant C . As a consequence, 

E{e LP ) < 00, 



for any p < -m, and there exists a constant Ao > such that E(e 
for all A < Aq. 



XL Hd 



) < OO 



Proof. From the equality 




and the definition of L e , we obtain 
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This expression allows us to compute the moments of L E . Fix an integer 
n > 1. Denote by T n the set {0 < s < t < T} n . Then 



E{L n £ ) 



jf jf ^ E [exp (i (^, Btl -B sl ) + ... + i (e„, fl f „ - B Sn ))\ 

(3.4) 

where s = (si, . . . , s n ) and t = (ti, . . . , t n ). Notice that 



exp ( ~2 ^ '^ ' 2 I d ^ " ' d £ ndsdt > 



/ £ [exp (i (a, 5 tl - S S1 ) + • - • + i (U, B tn - B Sn })] 

JR nd 



xe 2^ J= iM 



1 



d£x ■ ■ ■ di n 



exp --E 

nd \ Z 



xe 2 ^ 



^6 • • - rfCn 



(/ 



1 



(3.5) 



where Q is the covariance matrix of the n-dimensional random vector (.6/ 
B\ — B\ ). Substituting ([33]) into ([331) y ields 



(27T;— j Tr 

and E(L™) converges as e tends to zero to 

1 ttf ( 1 



a, 



(27T) 

1 



exp ( -\^Qt ) ^ ) dsdt 



(det Q) 2 dsdt, 

provided a n is finite. 

If «2 < °o, then in the same way as before we obtain 

lim E(L £ L S ) = a 2 , 

£,610 

which implies that L e converges in I? as e tends to zero. Furthermore, if a n 
is finite for all n > 1, then we deduce the convergence in LP for any p > 2 
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of L £ as e tends to zero. The limit, denoted by L, will be, by definition, the 
self-intersection local time of the process B in the time interval [0, T]. To 
complete the proof of the theorem it suffices to show that a n is bounded by 
C n (n\) Hd , for some constant C. 
We can write 

TI. f d 

a n = I (detQ)~2 dsdt. 

(2n)^~ JTuDih < -<t n } 

For each i = 1, . . . , n we denote by T{ the point in the set {sj, s»+i, . . . , s n , 
which is closer to ij from the left. Then, by (HI) and the fact that Sj < ti, 
i = 1, . . . ,n, we obtain, using Lemma [5] in the Appendix, 

detQ = VsxiBl-BljVaxiBl-BljBl-Blj 

x • • • x Var(^ n - BljBl - B x n ,. . . X_ a - 

> Var( J B t 1 i |£ s \ )Var(5j, |Bj , fl a \ , # a ) 

x ••• x Var « 1 4 , ^ , . . . , ^ , B] n _ x , < ) 

> Var^J^jVar^J^) • • • Var(^J F Tn ) 

> k?(h - r x ) 2H {t 2 - r 2 ) 2H ■■■(t n - r n ) 2H . 

As a consequence, 

a n < '-^k t 2 / \\{U - n)- Hd dsdt. 

(27r) 2 JT n n{ti<-<t„} i=1 

If we fix the points t\ < ■ ■ ■ < t n , there are 3 x 5 x • • • x (2n — 1) = (2n — 1)!! 
posible ways to place the points si,...,s n . In fact, si must be in (0,ti). 
For s 2 we have three choices: (0, s%), (si,ti) and (ii, *2)- By a recursive 
argument it is clear that we have (2i — 1) possible choices for Sj, given 
s\, . . . , In this way, up to a set of measure zero, we can decompose the 
set T n H {ti < • • • < t n } into the union of (2n — 1)!! disjoint subsets. The 
integral of Ui=i(ti ~ n)~ Hd on each one of these subset can be expressed as 

J {0<Zi<-<z 2n <T} i=1 

where <r(l) < • • • < cr(n) are n elements in {1, 2, ... , 2n}, and z = (zi, . . . , z 2n ) 
Making the change of variables Hi = Z{ — i = 1, . . . , 2n (with the con- 



S 



vention Zq = 0) we obtain 



** = Uy^dy^ Uyr Hd dy 

J{0<yi+-+y2n<T} i=1 n! i(0<sif4j„<T} j=1 

_ 1 T n(2~Hd)+Hd F(l - Hd) n ~ l 



Therefore 



rrtn 

n\ r(n(l- Hd)+Hd+1)' 



— — 

k x 2 (2n - i)!!T n ( 2 -^) +/M r(l - i^)™" 1 

a ™ — nd 

(2vr)-r(ra(l - ffd) + Hd + 1) 



CiC „. (2n-l)!! 



r(n(l - Hd) + Hd+iy 



with Ci = T Hd r(l - Hd)- 1 , and C 2 = 1 ( j } . Taking into 

(2tt)2 

account that (2n - 1)!! < 2 n ~ 1 n!, and that 

r(n(l - Hd) +Hd+1)> C n (n\) 1 - Hd , 

for some constant C, we obtain the desired estimate. ■ 

If Hd > 1, the above result is no longer true. In that case the expectation 
of L £ blows up as e tends to zero. In fact, if we denote a 2 (s,t) = Var(S t 1 — 
Bj), for s < t, then 

E{L £ )= [ [ p £+a 2 {s , t) {0)dsdt = (2tt)-3 / / {e + a 2 (s,t))-Usdt, 
Jo Jo Jo Jo 

which converges to 

rT ft , a f T ft 



(2tt)~5 f f a 2 {s,t)-idsdt>{2^)-ik l d2 [ f (t- S y Hd dsdt 
Jo Jo Jo Jo 



CO. 



In this case, one can study the existence of the renormalized self-intersection 
local time defined as the limit as e tends to zero of L £ — E(L £ ). In the next 
section we discuss the existence and exponential moments of the renormal- 
ized self-intersection local time, using Clark-Ocone formula, in the case of 
the fractional Brownian motion. 
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4 Renormalized self-intersection local time of the 
fBm 



The fractional Brownian motion on M. d with Hurst parameter H € (0, 1) is 
a (i-dimensional Gaussian process B H = {B?,t > 0} with zero mean and 
covariance function given by (jl.ip . We will assume that d > 2. 

It is well-known that B H possesses the following integral representation 



K H (t,s)dW s 



where W = {Wt,t > 0} is a d-dimensional Brownian motion, and Kjj(s,t) 
is the square integrable kernel given by 



K H (t,s) = C H ,is 7 - 



-H 



(U 



\H - 



2 U 



H- 



2 du. 



if H > i, and by 
K H (t,s) = C H ,2 



(t - sf~ 



Is 1 H 



(u 



s) H ~^du 



if H < 2, for any s < t, where the constants are Cjj. 



H(2H-1) 



and Ch ,2 



B(2-2H,H-±) _ 

, where B(a,f3) denotes th beta function. 



_(l-2H)b(l-2H,H+±) 

The processes B and generate the same filtration, that is, Tt = 
(t{W s ,0 < s<t} = a{B^,0 <s<t}. 

The fractional Brownian motion satisfies the following local nondeter- 
minism property: 

(LND) There exists a constant hi > 0, depending only on H and T , 
such that for any t G [0, T], < r < t A (T — t) and for i = 1, . . . , d, 



VarC^lflf'* 



\s — t\ > r) > k<i r 



2H 



(4.1) 



Consider the approximated self-intersection local time L e introduced in 
(|3.3p . From the general result proved in Section 2 it follows that if Hd < 1, 
then L £ converges in L 2 to the self-intersection local time L, and the random 
variable L has exponential moments. If Hd > 1, this result is no longer true, 
and one considers the renormalization of the self-intersection local time, 
introduced by Varadhan. 
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The purpose of this section is to apply the Clark-Ocone formula to pro- 
vide a stochastic integral representation for the renormalized self-intersection 
local time L. As a consequence, we will prove the existence of some expo- 
nential moments for the random variable L. 

Theorem 2 Suppose that H < min -^^j ■ Then the renormalized self- 
intersection local time of the A- dimensional fractional Brownian motion B H 
exists in I? and it has the following integral representation 

Z = -E/ / / ^P*l s Mr,t,s)l K H(t,r)-K H (s,r)]dsdt)dWi, 

~^ J0 \Jr JO V r ,s,t " J 

(4.2) 

where 

A r j, s = E(Bf - Bf\F r ) 

and 

al, = V a r(B^-B^\F r ). 
Proof. The proof will be done in several steps. 

Step 1 We are going to apply Clark-Ocone formula to the random vari- 
able L £ . It is clear that L £ belongs to B 1,2 , and its derivative can be com- 
puted as follows 

D\L e = £ jf ^-(Bf - Bf)I? r (B^ - Bf •*) dsdt, 

where r G [0, T], and i = 1, . . . , d. Using 

Dl {B^ - Bf'*) = [K H (t,r) - K H (s,r)] l [0)t] (r), 

we obtain 

DIL £ = £ J* ^-(Bf - Bf) [K H (t,r) - K H {s,r)\ dsdt. (4.3) 

The next step is to compute the conditional expectation E{D t r L £ \J : r ). The 
conditional law of Bf — Bf given T r is normal with mean A r ^, s and covari- 
ance matrix a^ st Id, where Id is the ci-dimensional identity matrix. Hence, 
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the conditional expectation E \7&:(Bf — B^)\!F T ) is given by 



, dxi 



f dp 



-(A 



o \-"-r,t,s) 
OXi 



r,s,t 

As a consequence, from (j4.3j) we obtain 

E {D*L s \Jy) = - f T f Ar '\ s Pe+a2 (A r , t , s ) [K H {t,r) - K H (s,r)] dsdt, 

Jr JO £ + &r,s,t 

and this leads to the following integral representation for L £ — E(L £ ) 
L £ ~ E{L £ ) 

~1 J0 \Jr Jo e + (7 r.s.t " I 



i=l 



Step 2 In order to pass to the limit as e tends to zero we proceed as 
follows. Set 



!]'(/-. /. .s) = f'\ s Pe+a 2 (A r , t)S ) r) - /■)] . (1. 1) 



Clearly, T, l e (r,t,s) converges pointwise as e tends to zero to 

A 



£ l (r,M) = ^ P(t2 (A r , t)S ) [^(t,r) -K H (s,r)] 



In order to establish the convergence of the integrals in the variables s and 
t, we will first decompose the interval [0, t] into the disjoint union of [r, t] 
and [0, r). In this way we obtain 

L e -E(L £ ) = L^+Lf\ 

where 

L e 1) = -fll [I J r K(r,t,s)dsdt)dW*, 
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and 

d pT / pT pr \ 
L z ] = ~ E J [J I S ^ r ' £ > S ^ dsdt ) dW r- 

Step 3 We claim that the random field S|(r, t, s) is uniformly bounded 
on the set 0<r<s<ibyan integrable function not depending on e. 
In fact, using the local nondeterminism property (LND), and Lemma [5] 
in the Appendix, we obtain the following lower bound for the conditional 
variance of s t = Var(£>^'* — B s ,i \J r r ): 

al s , t > Var(v3f ' l - B^\T S ) = Var(sf 'V.) > ~ s) 2H ■ (4.5) 

We can get rid off the factor A l rt s in the expression (|4.4|) of S* (r, t, s) using 
the inequality 

_d i 1 _d , 1 

t 2~ 2 \x\ z t 2" r 2 

Pt{x) < C——e-— < C——, (4.6) 
\x\ \x\ 

for some constant C > 0. In this way we obtain, using (j4.5j) and ([4,61) 

|S*(r, t, s)| < C (i - s)- Hd - H \K H (t, r) - K H {s : r)\ , (4.7) 
for some constant C > 0, and by Lemma [7] in the Appendix we obtain that 

/ / (t- sy Hd ~ H \K H {t,r) - K H (s,r)\dsdt <C(A~ H VI). (4.8) 

Jr Jr 

By dominated convergence we deduce the convergence of the integrals 

lim f f El(r,t,s)dsdt= f f T, i (r,t,s)dsdt 

for all (r, uj) £ [0, T] x $7, and a second application of the dominated conver- 
gence theorem yields that ff, Y, l £ (r,t, s)dsdt converges in L 2 ([0,T] x U) 
to f*'E i (r,t,s)dsdt. This implies the convergence of Le to 



Jr Jr 



I/ j r T ' i ^^' S ^ dsdt ) dW r 



in L 2 (f2) as e tends to zero. 

Step 4 Consider now the case s < r < t. In this case the integral 
of the term S|(r, i, s) is not necessarily bounded, and in order to show the 

(2) 

convergence of L e we will prove uniform bounds in e for the expectation 
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E §1 |Si(r, t, s)\ p dsdt^j , for some p > 1. We can write for s < r < t, 
using the first inequality in (|4.6p 



|S*(r,t,s)| < I^L Pe+< j (A r ,t, s )|^(t,r) 

I >T -4- <T . 1 r-,8,1 



< g^ + ^J^exp ^- 4( ^f t) ] \K H (t,r)\, (4.9) 

for some constant C > 0. If s < r < t, using the local nondeterminism 
property (LND) we obtain the following lower bound for the conditional 
variance a 2 s t : 

a 2 r s t = Var(Bf '* - B? '*|J>) = Var(sf '^Jv) > fc 2 (i - r) 2H . (4.10) 

On the other hand, if s < r < t 

a 2 ^ t = V a v(B?> i -B?> i \T r ) = Var(B?> i -B?' i \F r ) 

< Vax(B? '* - B? = (t - r) 2H . (4.11) 

Also we will make use of the estimate (see [8]) 

\K H (t, r)\ < k 3 (t - r) H -^A~ H . (4.12) 
Substituting the estimates KWh . (14TTTJ) and KT2b into (079]) yields 

|E*(r,*,s)| <CA- H ^ £ (r,t,s), (4.13) 
for some constant C, where 

tf e (r, i, a) = (e + & 2 (t - r) 2H ) (t - r) H ^ exp ( 1 ' ^ 



4(e + (i - r) 2 ") 

(4-14) 

Notice that if < i, then £*(r, t, s) is uniformly bounded by the 

integrable function Cr~2~ H (t — r)~ Hd ~^ , and we can conclude as in Step 
3. For this reason, we can assume that Hd > i. 
We claim that for some p > 1 , we have 

su P e([ [ (r, t, s)dsdt] < oo. (4.15) 

e>0 V^r JO / 



14 



To show this estimate we first derive a lower bound for the expectation of 

r i 2 

\Al : t,s\ 2 = E(B t ' — B s ' \T r ) ■ The main idea is to add and substract the 

term B^' 1 , and then neglect the expectation E ( (^E{B^ I,1 \ J- r ) — B,?' 1 ^ ^ ^ . 

This argument will be used later to find a lower bound for the covariance 
matrix of the vector (^E^B^' 1 — B^f' 1 ).? 7 ,,), 1 < i < nj . 



MI4mI 2 ) 



E ((e(B^ - B^\T r ))^ 
E ^E(B^\T r ) - B^ 1 ) 2 ^ 

+2E ((E(B^\F r ) - B?A (B^ 1 - Bf- 1 )) + E ( (B* 1 



> 



2E ((sf 1 - B?*) (B?' 1 - Bf + E ((Bf - 1 - B? 



L ( ( Bf - 1 - Bf - 1 

>2H 



B B 



(t - s) 2ii - (t - r) 2// . 
As a consequence, we obtain, assuming p < 2 



E ( exp ( — 



p|A-,*, s | 2 



4(e + (t - r) 2 ^) 



l + |( £ + (t-r)^)-iB(|A 



1 rY 



\2H 



Hence, 



B exp 



P|-4r,t, S | 2 



4(e + (t - r) 2 ^) 



< C(e+ (t-r)^)f(t-r)- 2 ^ a (t-s) 



-2/773 



(4.16) 
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where a + f3 = f . Substituting (f^7T6]) into P~T4"]) yields 



E (lr L q ' P ^ t '^ dsdt ^ ~ C J j Q ( £ + ( t - r Y 

x(i - r)( H -^) p (t - s)-P 2H dsdt 

< C [ T [\t- r )-PHd-i+2H0 {f _ s) -2H(3 dsdL 
Jr Jo 



If -ffci > 1, we can choose (3 such that 2ff/3 > 1, and integrating in the 
variable s, the above integral is bounded by 

C f T (t-r)- pHd -i +1 dt, 



which is finite it p > 1 satisfyes (Hd + ^) p < 2 (this is possible because 
Hd+ \ < 2). If Hd < 1, we can choose (3 such that 2H(3 = Hd — 5, for any 
5 > , and we obtain the bound 



■J r 



which is again finite if p > 1 is close to one, and 5 > is small enough. 

As a consequence, from (|4,13p and (|4.15p . for any fixed r £ [0,T], the 
family of functions { S*(r, i,s),e > 0}, is uniformly integrable in [r,T] x 
[0, r], so it converges in L l ([r, T] x [0, r]) x f2 to X*(r, i, s), for i = 1, . . . , d. 
This implies the convergence of the integrals 

pT pr pT pr 

lim / / T,i(r,t,s)dsdt = / / E*(r, t, s)dsdt, 

eiO Jr Jo Jr Jo 

for each fixed r £ [0, T] in 

Finally, we claim that this convergence also holds in L 2 ([0,T]xO), and 

(2) 

this implies the convergence of L e to 



in L 2 (Q) as e tends to zero. To show the convergence in L 2 ([0, T] x Q) of 
the integrals 



j-T pr 

Y*(r) = / T,i(r,t,s)dsdt 

Jr Jo 
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it suffices to prove that 



sup / E (\Y £ i (r)\ P ) dr < oo 

£>0 Jo ^ ' 



(4.17) 



for all i = 1, . . . , d and for some p > 2. The proof of (|4.17p will be the last 
step in the proof of this theorem. 

Step 5 Suppose first that Hd < 1. Then, from (|4. 13[) we obtain 



i-T . r T 



E |F e l (r)r dr<C E 



t-T t-r 
K JO 



V 

ty £ (r,t,s)dsdt I 



Using (j4.14H and Minkowski's inequality yields 



ty £ (r,t,s)dsdt 







< 


ir So 


V 





£ + k 2 {t-r 



<2H 



cxp 



\Ar, 



t.s 



and from (|4. 16[) . choosing (3 = &, we get 



cxp 



\A + I 2 



4(e + {t- rf H ) 



<C{e+(t-r) 2H )^(t-s)~T. (4.19) 



(t-r) H ~2 
dsdt, (4.18) 



4(e + (t - r) 2H 

Substituting (|4TT9|) into (I4TT81) yields 
/•T /-r 
Ir Jo 



^ e (r,t,s)dsdt 



< C [ T (t-r)' Hd ~^ +I Tdr, 



which is finite if we choose p > 2 such that p < 2 Hd~i • Fi nanv > if P (3 ~~ -H") > 
— 1 we complete the proof of (|4.17|) in the case Hd < 1. 

In the case iifd > 1 we cannot apply the previous arguments, and the 
proof of (|4.17p follows from the moment estimates given in Proposition [3j 



Remark 1 Theorem [2] also provides an alternative proof of the existence of 
the self-intersection local time in the case H £ [^,min(^, which was 

proved by Hu and Nualart in [9] in the general case Hd < |. Notice that for 
d > 3, the condition H € [1, min(^, ^-)) is equivalent to 1 < Hd < |, and 
for d = 2 we require H < instead of the more general condition H < |, 
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that guarantees the existence of the renormalized local time (see |16j and 

The next Proposition contains the basic estimates on the moments of the 
quadratic variation of the stochastic integral appearing in the representation 
of the renormalized self-intersection local time. 

Proposition 3 Assume l<Hd<\. Set 



A £ (r) 



T r r 



r JO 



<& £ (r,t,s)dsdt, 



where *$> e (r,t,s) has been defined in {^.1$ . Then, for any integer n > 1, 

E(K n £ (r))<C n (n^, 
for some constant C > 0, where 



7 > 



1 



+ H\\d 



1 
2H 



Proof. Set g £ (t - r) = [e + k 2 (t - r) 



2H 



(t - r) 2 . We have 



E{K(r)) = E 



I Jo 9e<(t ~ r ' >eXP ( 



14 J 2 



4(e + (t - rf H ) 



dsdt 



r,T]n J Sn f = \ 



\A l I 2 



^(e + iU-r)™) 



dsdt, (4.20) 



where S n = {0 < si < ■ ■ ■ < s n < r}, s = (si, . . . , s n ) and t = (ti, . . . , t n ) 
We denote by Q the covariance matrix of the vector 

{e(B*> X - B**\F r ), E(B% 1 - Bf> Vr)) • 
Then, a well-known formula for Gaussian random variables implies that 



E 



exp ( - 



\A l I 2 



^ 4(e + (U - rf H ) J 



det( I + -QD- 1 



2^ fl^El det (2D + Q)~^ , (4.21) 
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where D denotes the n x n diagonal matrix with entries Oj = e+ (U — r) 2H . 
As in the computation of E(jAl ts \ 2 ), adding and substracting the term 



B^ 1 ' 1 yields 



Qij = E(E(B^ - B^ l \T r )E(B^ 1 -B^\T r 



= E (E(B^ - B^ l \F r )E(B^ - B^ x \T r 

+E ([B?* - B^)(B^ - B?*)) + E ((Bf 1 - B^)(B^ - B^)) 

+e [(B^ - B Si )(B^ - B^)) 
= E (E(BP - B^\F r )E{B?; 1 - B^\T r )) 

-E ((<• 1 - B^XB*' 1 - B^)) + E ((Sf 1 - B^)(B^ - B?; 1 )) . 

Hence, we obtain 



Q = R- N + M, 



where 



Rij = E{E{B^-B^ r )E{B^-B^\T r )), 
Ma = E ((Bf 1 - B^)(B^ - Bj 1 )) , 
N t3 = E ({B*' 1 - B?' 1 )^' 1 - B?' 1 )) . 

All these matrices are nonnegative definite. The main idea will be to get 
rid off the matrix R, and control the matrix N by its diagonal elements 
which are 

N u = (U-r) 2H . 

Indeed, the matrix N is nonnegative definite and, hence, it safisties the 
inequality 

N < nD N , (4.22) 
where Djy is a diagonal matrix whose entries are Na. Therefore, 

Q > -N + M > -nD N + M, 

and for any 1 < 6 < 2, we can write 

det(2D + Q) > det(2D + ^—^-Q) < det(2D-(2-S)D N + ^—^-M). (4.23) 

n n 
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The entries of the diagonal matrix D\ = 2D — (2 — 5)Dn are the positive 
numbers 

2e + S(U - r) 2H > 0. 
Prom (jOOj) . (|43T|) and P~23j) we obtain 

£(A"(r)) < 2^n! / / f[ (g E {ti ~ r)af) 

J[r,T] n JS n V / 

2 S d 

x det(L>i + M)~2dsdt. 
n 

We have 

det(Di + ^—M)-i < ( ) (det Z?i)" Q (det M)~ p , 

n \2 — 6 J 

where a + f3 = 77. Hence, 

2^nl / ]jU(t i -r)af(2s + S(t i -rf H y a ) 

X (det M) _/3 dscft. 

Then, 

g £ (ti-r)aj (2e + 2(t { - r) 2H y a 
= (e + hiU-r^y 2 {tl -r) H ^{e + (t l -rf H y{2e + 2(t l -rf H )- a 

< C(ti-r)-^- 2Ha , 
for some constant C > 0. Thus 

E(A£(r)) < C n n^ n n\ f [ Y\(U - r)~^~ 2Ha (det M^^dsdt, (4.24) 

J[r,T] n Js n fj[ 

for some constant C > 0. 

Applying Lemma[5]in the Appendix and the local nondeterminism prop- 
erty of the fractional Brownian motion we obtain 

detM = Var(£ tn - B Sn )Var(B t „_ 1 - B Sn _,\B tn - B S J 
x • • • x Vax(B tl - B Sl \B t2 - B S2 , . . . , B tn - Bj 

= (t n ~ 8 n )* B Vax(B an _ 1 \B tn _ l ,B tn ,B an ) 
x ••• x Vax(B sl \B tl ,...,B tn ,B sl ,...,B Sn _ 1 ) 

>kr\r - s n ) 2H {{s n - Sn _0 A s^) 211 ■ ■ ■ ((s 2 - Sl ) A Sl ) 2H . (4.25) 
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Substituting (|4.25j) into (|4.24f) . and choosing a such that a < -^j (this is 
possible because Hd > 1) yields 

E (A e ») < C n n^n\ f \{r - s n ) ((s n - s n ^) A s^) ■ ■ ■ ((s 2 - si) A si)]" 2 ^ ds. 



Finally, by Lemma [5] in the Appendix we obtain 

C n n^ n n\ 



r(n(l - 2HP) + 1) 
Notice that /3 = | — a > | — And hence, 

E(A£(r)) < C n (n\ p + 2Hf *, 

where 

W+ 2H)>|- J L +M -i=(i + ff )( d -^ 
This concludes the proof. ■ 

Using the above proposition we can deduce the following integrability 
results for the renormalized self-intersection local time. 

Theorem 4 Assume ^ < H < min -^^j ■ For any integer p < ^ [(^ + if) (d — t^)] 1 
we have 

£(exp|L| p ) < oo. 

Proof. Taking into account Lemma [6] in the Appendix, it suffices to 
show that 

E ( exp ( L) ) < oo, 



where 

^) = Ylf (/ J^^ i (r,t,s)dsdtj dr. 
As in the proof of Theorem [2] we make the decomposition 

\ T l i (r,t,s)dsdt = / / E*(r, t, s)dsdt + / / H l (r,t, s)dsdt. 

r JO Jr Jr Jr Jo 

From (|4T7|) and (|4.8p we know that 



Z l (r,t,s)dsdt 



< C{A~ H V 1). 
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Therefore, applying Fatou's lemma and the estimate (|4.13p yields 



£(exp(L) ) < CE\ exp 



d „t 




E 

i=l 



'/' /■/' \ 2 

S*(r, i, s)dsdt J dr 



< C lim inf E exp 



T JO 



E 



i=i 



< C lim inf £ I exp I C 



o 



T 



T ,-r \ 2 

Z,i(r,t,s)dsdt ) dr 



r JO 



A-2H 







Applying Holder and Jensen inequalities we obtain 



£(exp ( L 



< C lim inf E 1 I exp I C 



A-2H 







< C lim inf 

^40 J 



,l-2H 



E exp C 



T 

r JO 



^ £ (r,t,s)dsdtj dr 



T ,-r \ 2 'P N 

^ £ (r,t, s)dsdt ) cir 



r JO 

T 



^> e (r,t,s)dsdt 



2p\ 



dr. 



Finally, 



E exp C 



n=l 



T rr 



2p\ 



r JO 
T 



^ £ (r,t,s)dsdt I 



r JO 



*$> e (r,t,s)dsdt 



2np\ 



oo „ n 

< —mn P ] + r 

L / 77 



ivy 



n=l 



and it suffices to apply Proposition [3] to conclude the proof. ■ 

Remark 2 The exponent po = \ \i\ + H) ~~ Jh)] * ^ s no * optimal. For 
instance, if Hd = 1, then p = and we know that for Hd < 1, then 

P° = TJd- ^ n particular, if H = ^ and d = 2 we obtain po = \i an d we know 
that in this case the critical exponent is po = 1- The lack of optimality 
is due to the factor n in the estimation of the positive definite matrix N 
by its diagonal elements given in (|4.22p . Without this factor n we would 
get the critical exponent 2 Hd-\ ' Du ^ our me thod does not allow to get this 
value. 
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Remark 3 In the case of the planar Brownian motion B = {Bt,t > 0} 
(that is, d = 2, and H = formula (|4.2|) yields 



(4.26) 



The quadratic variation of this stochastic integral is 



2 



£ S? f (f I T~?F exp i J -W^f) isit ) * 



exp I — — ) ds ) dr 



vr 2 7o Wo \B r ~B s \ 2(r-r 



(X r — b r ) , 



71-2 Jo \Jo \B r — B s 
From Ito's calculus we know that 

rr ds 1 



/ o I B r — B s | d — 1 

where X r has the law of the modulus of a d-dimensional Brownian motion 
at time r (Bessel process), and b r has a normal N(0,r) law. We can write 

dr, 




which clearly imply the existence of some Ao such that E ^exp ^A (^Ljjj < 
oo for all A < Ao- From Lemma Owe get that there exists (5q such that 
E ( exp f/3 



L J J < oo for all /3 < /3q. This method does not allows us to 
obtain the critical exponent, just the existence of exponential moments. 

Remark 4 The above results remain true if we replace the fractional 
Brownian motion with Hurst paramter H, by an arbitrary centered Gaus- 
sian process of the form (|3.ip satisfying the local nondeterminism property 
(LND) and following properties: 

(CI) For any s,t £ [0, T], s < t, there exist constants k% and fc 4 such that 
k 3 (t - sf H < E{\B\ - B\\ 2 ) < k 4 {t - s) 2H . 
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(C2) The kernel K(t, s) satisfies the estimates 

\K(t,8)\<k 5 (t-8) H -*8*- H , 

for all s < t, and 

f [ (t- S y Hd - H \K{t,r)-K{s,r)\dsdt< ip(r), 

J r J r 

where J Q T ip(r) 2 dr < oo. 

5 Appendix 

In this Appendix we will first state and prove some elementary lemmas. The 
first one is well-known. 

Lemma 5 Suppose that Q\ C Q2 are two a-fields contained in T '. Then, for 
any square integrable random variable F we have 

Var(F|&) > Var(F|g 2 ). 

Let M = {M t , t > 0} be a continuous local martingale such that Mq = 0. 
Then, the following maximal exponential inequality is well-known 

P ( sup \M t \ > S, {M) T <p) < 2exp (~ 

\0<t<T J V 2 P 

As a consequence of this inequality we can obtain exponential moments for 
Mt from exponential moments of the quadratic variation (M)t 

Lemma 6 Suppose that for some a > and p G (0, 1] we have E(e a ( M ^T ) < 
00 . Then, 

° A I Mt I "> 



(i) if p = 1, for any A < y — , E(e ' T ') < 00, and 
(ii) ifp < I, E{e x \ MT \ P ) < 00 for all A > 0. 
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Proof. Set X = \Mt\ p ■ For any constant c > we can write 

/>oo 

E(e xx ) = / P(X > y)\e Xy dy 
Jo 

= / [P(X > y, (Mf T < cy) + P(X > y, (M) P T > cy)] Xe^dy 
Jo 

< ^°°2ex P ^-^j ^ Xvd v + l°° p (^— > y) XeXyd y 

= [°° 2A exp (\y - ^] dy + E(ec^T ) . 
Jo \ 2cp J 

Then it suffices to choose c = ^ to complete the proof. ■ 

The next two results are technical lemmas used in the paper. 

Lemma 7 Suppose that H < min(^-j-, ^). Then, we have 

£ J\t - sY Hd - H \K H (t,r) - K H (s,r)\dsdt < C (r*~ H V l) , 

for some constant C . 

Proof. We know that 



2 , 



Then 



/:= f f (t - s)- Hd - H \K H (t,r) - K H (s,r)\dsdt 

J r Jr 

< c£ J* (t - s)- Hd - H (?) H ~ (9 - rf-ldOdsdt. 



If H < i, then, » < 1, and if H > \, then (^) H * < cA~ H . Hence, 

the above integral is bounded by 



C{A~ H V 1) f T f f (t - s y Hd - H (9 - r) H -ld9dsdt. 

Jr Jr J s 
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From the decomposition 



~-H = a + (3, 
Hd + H = 7 + 5, 



we obtain 

i-T rt rt 



r J r J s 
T r-t r t 



r J r J s 



(t - s y Hd - H {9 - r) H -2d9dsdt 

(s - ry a {9 - s)-P-~<(t - 9)- & d9dsdt. 



Finally, it suffices to show the parameters a, /3, 7 and 5 in such a way that 
a < 1, 5 < 1 and (3 + 7 < 1. This leads to the condition 

1 3 

-+Hd< min(l, - - H) + min(l, Hd + H), 
which is satisfied if H < min(^-, ^). ■ 

Lemma 8 Let a < 1. Fix an interval [0, T]. For each integer n > 1 we 
have 

[((T - s„) A s n ) ((s n - s n _i) A s n _i) • • • ((s 2 - si) A si)p a (is 

A n (T) 

rpn(l— a) 

< x ;C n , (5.1) 

-r(n(l-a) + l) V ; 

where A n (T) = {0 < si < ■ ■ ■ < s n < T} 

Proof. We proceed by induction on n. For n = 1 we can write 

((T-sx) A Sl y a d Sl = 1 2 s^ a d Sl + [ (T- Sl )- a d Sl 



2 /T x 



1 - a \ 2 

which implies (J5H]) with C = ^^-2 a . 
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Suppose that the result holds for n — 1. Then, 

I n = [((T-s n ) As„)((s„-s n _i) A s„_i) ■ ■ ■ ((j»2 - si) Asi)] _a ds 

JA n (T) 

= / ((T-s n ) As n )-° 
Jo 

x / [((s„ - s n _i) A s n _i) • • • ((s 2 - si) A si)] _a dsi ■ ■ • ds n -i 

\J A n _i(s n ) 

By the induction hypothesis we can write 

fin— 1 /»T 

/n < -J ((T- Sn )As n )- a sir m - a) ds n 

T(n - a) J 



C 



n—l 



r((n-l)(l-o) + l) 

x I J 2 s^- 1 )( 1 -°)- a ds n + j T {T-s n )~ a s^' l ^ l - a Us n 



< 



C 



2 

n—l 



T(n(l - a) + a) 



n(l - a) \2 



n(l— a) 



^ 1 (l-^)" a ^ (n_1)(1_a) ^ 



< r^-^c- 1 / i r(i-a)r((n-i)(i-o) + i) 



T(n(l - a) + a) \n(l - a) T(n(l - a) + 1) 
^(1-0)^-1 f I , r (! ~ Q ) 



n{l - a)r(n(l - a) + a) T(n(l-a) + l), 
Using the relation r(n + 1) = nT{n) we obtain 

rc(l - a)r(n(l - a) + a) > n(l - a)r(n(l - a)) = T(n(l - a) + 1), 
and, as a consequence 

I n < T^C"- 1 (1 + T(l - a)) 1 

T(n(l - a) + 1) 

and it suffices to take C > max {^^-2 a , 1 + T(l - a)) . ■ 
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